HURWITZIAN CONTINUED FRACTIONS CONTAINING A 
REPEATED CONSTANT AND AN ARITHMETIC PROGRESSION 



GABOR HETYEI 

Abstract. We prove an explicit formula for infinitely many convergents of Hur- 
witzian continued fractions that repeat several copies of the same constant and ele- 
ments of one arithmetic progression, in a quasi-periodic fashion. The proof involves 
combinatorics and formal Laurent series. Using very little analysis we can express 
their limits in terms of (modified) Bessel functions and Fibonacci polynomials. The 
limit formula is a generalization of Lehmer's theorem that implies the continuous frac- 
tion expansions of e and tan(l), and it can also be derived from Lehmer's work using 
Fibonacci polynomial identities. We completely characterize those implementations 
of our limit formula for which the parameter of each Bessel function is the half of an 
odd integer, allowing them to be replaced with elementary functions. 



Introduction 

It is a remarkable property of infinite continued fractions that they often define a 
sequence whose limit is easier to describe than the individual entries. Even for [1,1,...], 
the simplest example, the limit is easily found by solving a quadratic equation, it takes 
a bit longer to find a formula for the convergents, in terms of Fibonacci numbers. 

The class of Hurwitzian continued fractions of the form 



a, /3 , a, . . . , a, f3 + & x -n 
^ .. ' 



oo 

71=1 



that we propose investigating seems to be no different in this regard. For the special case 
r = 0, D. H. Lehmer [S] proved a formula for the limit in terms of Bessel functions, that 
can be verified very easily, after having guessed the correct answer. Using another result 
from [9] and some Fibonacci polynomial identities, it is not hard to generalize Lehmer's 
formula to Hurwitzian continued fractions of the above form (see Remark I2.10p . The 
resulting generalization has several famous special cases, the most famous ones being 
Euler's formula for e and the formula for tan(l). On the other hand, the only research 
regarding the convergents themselves seems to be the work [S] of D. N. Lehmer (D. 
H. Lehmer's father!), who proved congruences for their numerators and denominators 
by induction, in somewhat greater generality, but without stating the values of the 
convergents in an explicit fashion. 

Our work contains such an explicit formula, stated in Section [2] and proved in Sec- 
tion HJ A variant of the Euler-Mindig formulas using shifted partial denominators yields 
a summation formula with many vanishing terms. This leads to a compact recurrence 
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for the numerators of the (nd + r — l)st convergents that may be restated as a lin- 
ear differential equation for a formal Laurent series. Classical textbooks on differential 
equations instruct us to solve the associated homogeneous equation first and then find 
the general solution by "variation of parameters" . Unfortunately, in our case the "so- 
lution" of the homogeneous equation turns out to be a a two-way infinite formal sum. 
After discarding infinitely many terms to have a formal Laurent series, the "spirit" of 
the classical method still inspires a good guess for the form of the solution, where the 
transformed differential equation encodes a recurrence that is easily solved by inspec- 
tion. Returning to the original Laurent series involves using two polynomial summation 
formulas that can be shown purely combinatorially, and seem to be interesting by their 
own right. These formulas are shown in Section |3j As outlined in Section EJ our explicit 
formulas for the (nd + r — l)st convergents allow us to calculate the limits using very 
little analysis. To find only the limits, this approach is a bit more tedious than the 
one outlined in Lehmer's work [9], but we gain a little more insight by also obtaining 
asymptotic formulas for the numerators and the denominators of the convergents. 

Finally, Section [5] is a reflection on Komatsu's recent remark [7] stating that all 
known examples of Hurwitzian continued fractions seem to have a short quasi-period 
and involve (hyperbolic) trigonometric functions. For our class of Hurwitzian continued 
fractions, the limit may be expressed in terms of (modified) Bessel functions, which are 
known to have an elementary form, in terms of rational and (hyperbolic) trigonometric 
functions, when their parameter is the half of an odd integer. We describe all con- 
tinued fractions in our class that yield such parameters, and find that, subject to this 
restriction, the quasi-period d can not be longer than 4. 

Our work inspires several questions worth exploring in the future. Extending the 
validity of our formulas to all convergents of the same class seems to require only a 
little more work. As already indicated in Lehmer's paper [9], the calculation of the 
limits is easily extended to the class of Hurwitzian fractions whose quasi-periodic part 
contains several different constants and one arithmetic sequence. An explicit formula for 
the convergents should be obtained using some multivariate generalization of Fibonacci 
polynomials. The ideas used in our calculations of the convergents may also be useful 
in finding convergents of other Hurwitzian continued fractions. Finally, as indicated in 
Sections [3] and IH there is combinatorics behind the formulas for the convergents. This 
combinatorics would be worth uncovering. A good starting point may be revisiting 
the weighted lattice-path model proposed in Flajolet's work [3], where convergents of 
generalized continued fractions arise as weights of infinite lattice paths of limited height. 

1. Preliminaries 

1.1. Simple continued fractions. A good reference on the basic facts of the subject is 
Perron's classic work [13] [H]. A (generalized) finite continued fraction is an expression 
of the form 

h 

a + 7 , (1-1) 

ai H 

o n -i 

a 2 + ■■ — 

a n -i H 
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where the initial term do, the partial denominators and the partial numerators 

b\, . . . , b n may be numbers or functions. An infinite continued fraction is obtained by 
letting n go to infinity The arising questions of convergence have a reassuring answer 
for (simple) continued fractions where all hi equal 1, the initial term do is an integer, 
and the partial denominators Oj (for i > 0) are positive integers. In the present work 
we will only use simple continued fractions, and refer to them as continued fractions. 
Every rational number may be written uniquely as a finite continued fraction, subject 
to the restriction that the last partial denominator is at least 2 [T3| §9, Satz 2.1]. 
Every infinite continued fraction converges to an irrational number and every irrational 
number may be uniquely written as a (necessarily infinite) continued fraction [T3"l §12, 
Satz 2.6]. As usual, for finite, respectively infinite, continued fraction we will use the 
shorthand notations [do, a±, . . . , a n ] and [do, di, . . •], respectively. An infinite continued 
fraction [do, di, . . .] is the limit of its convergents, that is, of the finite continued fractions 
[do, . . . ,a n ], obtained by reading the first n partial denominators. The convergents 

[d , . . . ,a n ] = — (1.2) 

Qn 

may be recursively computed from the initial conditions p_\ = 1, p = a , q_i = and 
q = 1, and from the recurrences p n = a n p n _i + p n -2 and q n = a n q n _\ + g n _ 2 for n > 1, 
cf. [131 §2, Eq. (12), (13)]. The integers p n and q n are relative prime for all n [T31 §9, 
Satz 2.1]. The Euler-Mindig formulas, derived for general continued fractions in [T3(, 
§3], allow to express the numerators p n and the denominators q n directly. 

Definition 1.1. We call a set S of integers even if it is the disjoint union of intervals 
of even cardinality. Given two sets of integers S and T such that S C T , we say that 
T evenly contains S , denoted by S C e T or T Z> e S , if T \ S is an even set. 

Thus, for example {5} C e {1, . . . , 7} since the difference {1, 2, 3, 4, 6, 7} is the disjoint 
union of {1, ... , 4} and {6, 7}, whereas {4} $Z e {1, . . . , 4} as {1, 2, 3} can not be written 
as the disjoint union of intervals of even cardinality. The Euler-Mindig formulas for 
(simple) continued fractions may be restated as 

Pn = \ ai and qn= ( L3 ) 

5C e {0,...,n} ieS SC R {l,...,n} ieS 

The following observation is an immediate consequence of Eq. ( II. 3p . 

Lemma 1.2. The denominator q n associated to [do, di, . . .] is the same as the numerator 
Pn-i associated to [d 1; a 2 , ■ ■ ■]■ 

1.2. Hurwitzian continued fractions. A Hurwitzian continued fraction is a con- 
tinued fraction of the form [do, . . . , d^, 4>o(X), ■ ■ ■ , 0fc_i(A)]^ =o , where <f> , . . . , <fik-i are 
polynomial functions that send positive integers into positive integers. The definition 
given in [T3| §32] is easily seen to be equivalent. In the most trivial examples of a 
Hurwitzian continued fraction all functions <pj are constants, we then obtain a periodic 
continued fraction. The set of real numbers represented by a periodic continued fraction 
is exactly the set of quadratic irrationals [131 §20]. It should be noted that quadratic 
irrationals are sometimes excluded from the definition of Hurwitzian continued fractions 
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in some recent papers [7j, by requiring that at least one of the repeatedly used polyno- 
mial functions be non-constant. A real number is Hurwitzian if its continued fraction 
representation is a Hurwitzian continued fraction. A famous Hurwitzian number is 

e= [2, 1,2, 1,1, 4, 1,1, 6, 1,1, 8, 1,1, 10,...], 

see (131 §34, Eq. (10)]. The main result on Hurwitzian numbers is Hurwitz theorem [131 
§33] stating that for a Hurwitzian number £ £ R, and any rational numbers a, b, c, d, 
satisfying ad — be ^ 0, the number (a£ + b)/(b£ + d) is also a Hurwitzian number. 
Hurwitz theorem provides also some estimate on the degrees of the polynomial functions 
appearing in the continued fraction representation of (a£o + b)/ (fc£o + d), remains silent 
however on the issue how the length of the period may be affected by the transformation 
£o l— (a£o + b)/(b£o + d). A generalization of Hurwitz theorem may be found in [T6j . 

Hurwitzian numbers may be computed from their continued fraction representation in 
some special cases, the expansion of the set of examples is subject of ongoing research. 
A frequently overlooked first attempt may be found in Perron's book [T4"l §48, Satz 
6.3] which states (without proof) a formula for generalized continued fractions of the 
form (jl.ip having the property that the partial numerators b\ are all equal, and that 
the numbers a; form an arithmetic sequence. For simple continued fractions, Perron's 
result gives 

oo 

[ft, ft + ft, ft + 2ft, ft + 3ft, . . .] = ft Z U ■ (1-4) 

i 

to /^n!r(^+n+l) 

The same class of Hurwitzian continued fractions was revisited by D. H. Lehmer [9], 
who proved the following formula. 

[ft, ft + ft, ft + 2ft, ft + 3ft, . . .] = 7 /3o//3l _ 1 (2ft- 1 )/^ o//3l (2ft- 1 ) (1.5) 

where I v {z) is the modified Bessel function of the first kind 

4(2) = E r ( m+1 ) r(z/ + m + 1 )- ( L6 ) 

m=0 ' 

D. H. Lehmer's result was generalized in recent papers to other Hurwitzian continued 
fractions [51 [71 [11] . Some analogous results were found for Tasoev continued fractions 

El HH [12], defined as continued fractions of the form [a ; a k , . . . , a h J ^ =1 . As pointed 

m 

out by T. Komatsu [7J, most of these recent results involve Hurwitzian and Tasoev 
continued fractions where the length of the quasi-period does not exceed 4, Komatsu's 
work [TJ contains some sophisticated examples of longer quasi-periods. 

1.3. Bessel functions. Besides the modified Bessel functions I v (x), our formulas will 
also involve the (original) Bessel functions of the first kind 

Mz) ^ o r(m + i)r(v + m + i)- {L7) 
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There is an elementary expression for the functions I v (z) and J u {z) respectively, when- 
ever v is the half of an odd integer. Indeed, using T(l/2) = yfn and Y{z + 1) = zT(z), 
it is easy to derive directly from the definitions (11.61) and 1 1 1.7ft that we have 



J-cosh(z), hn(z) = \— smh(z) 

V 7TZ \ 71Z 




-1/2{Z) = \ I [l.l 

[2 

J_i/ 2 u) = \ — cos (2;), and Ji/ 2 u) = \ — sin(z). (1.9) 
V ttz \ nz 

(see pUl List of formulae: 44, 48, 182, 186]). The functions I v (z) and J v (z) satisfy very 

similar recurrence formulas: 

= I v -i{x) - —I v {x) and J u+1 (x) = -J u ^(x) + —J u (x), 

X X 

see [TJ 9.1.27, 9.6.26]. These allow us to find explicit elementary expressions for I u (z) 
and J u (z), whenever v is the half of an odd integer. In particular, for v = 3/2 we obtain 

/~2~/ sinh(x)\ ^2~ /sin(x) , A 

4/2 1 = \ — cosh(x) and J3MX) = \ — cos(x) , 

\ 7TX \ X J \ TTX \ X J 

(1.10) 

see pUl List of formulae: 45, 183]. 

1.4. Fibonacci and Lucas polynomials. The Fibonacci and Lucas polynomials are 
g-analogues of the usual Fibonacci and Lucas numbers. 

Definition 1.3. We define the Fibonacci polynomials F n (q) and the Lucas polynomials 

L n (q) by the initial conditions F (q) = 0, Fi(q) = 1, L (q) = 2 and L\(q) = q; and 
by the common recurrence P n {q) = q • Pn-i(o) + Pn-2(o) for n > 2, where the letter P 
should be replaced by either F or L throughout the defining recurrence. 

Fibonacci and Lucas polynomials are widely studied, they even have their own Math- 
world and Wikipedia entries. As a sample reference on Fibonacci polynomials see Yuang 
and Zhang [17], for further generalizations, see Cigler [2]. 

The coefficients of the Fibonacci and Lucas polynomials are listed as sequences 
A102426 and A034807 in the Online Encyclopedia of Integer Sequences [15]. In a 
recent work of Foata and Han [4] the coefficient of g n + 1_m in F n+ i(q) appears as the 
number of t- compositions of n into m parts. Here we record another combinatorial 
interpretation. The coefficient of g n + 1_m in F n+ i(q) is the number of m-element even 
sets contained in {1, . . . , n}: 

0C e SC{l,...,n} SC e {l,...,n} 

Equation ( II. lip may be shown by induction on n, using the defining recurrence. Just 
like for the usual Fibonacci and Lucas numbers, a closed formula for F n (q) and L n (q) 
may be obtained after solving the characteristic equation. We have 

F n (q) = _l=(pf-^) and (1.12) 
vV + 4 

L n (q) = P1+P2 (1-13) 
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where 



q + yV + 4 g - yV + 4 
Pi = ^ and p 2 = ^ . (1.14) 

2. Our main result 
All our results will be about Hurwitzian continued fractions of the form 



£(a, /?o, AL,d,r) := [ a, . . . , a, /3 , a, . . . , a, /3 + & ■ raj 



oo 
„ n=l 

</ 



where a, (3q, Pi and d are positive integers and and r is a nonnegative integer. In order 
to state them we need to introduce two magic numbers associated to such a continued 
fraction. 

Definition 2.1. The magic sum associated to /3q, d, r) is the sum 

o\pt,Pa,P\,d,r) = — h 7r~ET? \ i 

Pi P\Fd{ot) 

the magic quotient associated to £(a, /? , Pi, d, r) £/ie quotient 



p(<*,Po,Pi,d,r) 



Pi ■ F d (a) 



Whenever this does not lead to confusion, we will omit the parameters (a, /3q, Pi, d, r) 
and denote the magic numbers simply by a and p, respectively. Note that a does not 
depend on r and that p depends only on a, /3i and d. Our main result is the following. 

Theorem 2.2. The {nd + r — l)s£ convergent of £(a, /3 , Pi,d,r) may be written as 
p nd+r _i/q nd+r _i where p nd+r ._i flwen fry 

Pnd+r-i „ , (n - fc)! /n + o- - 1 - k\ k d _ r 

F r+1 (a) ^ ,, _ 9 , P + (-1) ^-r-i(a)- 



F d (a)«3? y ^ fc! V n-2k 

L(n-1)/2J 



in — fc — ljl/n + cr — 1 — fc \ 



fc! V n -2k- 1 

fc=0 v 



fc+i 



and q n d+ r ~i is given by 



q.nd+r-i „ f v ^ (n - fc)! /n + a - 1 - fc^ k > w _ r , , 

L(n-l)/2j i 



„. . „ (n — fc — l)!/n + (T — 1 — k\ 

fc=0 v 7 



p 



We postpone the proof of Theorem 12.21 till Section HI In this section we only show 
how to find an asymptotic formula for the numerator and the denominator, which yields 
essentially the same limit formula that can also be derived from Lehmer's work [9]. 
This is the only part of our paper where the notion of limits and convergence from 
analysis will be used, the proof of Theorem 12.21 will only involve purely algebraic and 
combinatorial manipulations. We begin with estimating the magic sum a. 
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Lemma 2.3. The magic sum a is always positive. 

Proof. Since (3q and (3± are positive, it suffices to prove that Ld(a) > aFa(a) holds for 
all positive integer a and all nonnegative integer d. This may be shown by induction 
on d, using the common recurrence of the Fibonacci and Lucas polynomials. □ 



Corollary 2.4. The falling factorial (a - 



■n■ 



(a+n— 1) • • • (a) is always positive. 



As a consequence, we may divide p n d+ r -i/ (Fd(a) n (3i), as well as qw+r-i/ '(Fd(a) n f3™) 
by (cr + n — l) n and ask whether these quotients converge as n goes to infinity. It will 
turn out that they do. 

Lemma 2.5. Let a be any positive real number and p be any real number. Then 

\n/2\ 



lim - — 

n— >oo (u 



1 



n 



^ (n - k)\ 

k=0 



k\ 



n 



-a-l- 

n-2k 



P 



E 

m=0 



m!(cr + m — 1), 



Proof. Note first that the series on the right hand side is absolute convergent for any 
p as each of its term has smaller absolute value than the corresponding term of e p ^ = 
J2m=o P m °~~ m l m ^-i an absolute convergent series. For any e > there is an M such that 



E 



m\(o~ + m 



< 



(2.1) 



m=M+l 

Let us consider now the expression of n on the left hand side. It may be rewritten as 



\n/2\ 

E 



P 



k=0 



k\(a + k-l), 



■ d 



where 5. 



in 



n,k 



k)i 



o- 



n 



The factor 8 n ^ may be estimated as follows: 

k 



n 



2k 



a 



n 



k 



in 



k) 



in 



2k +\) 



< 



n 



k 



a + n — 1 



{a + n — 1) • • • {a + n — k) 

Our upper bound for 5 n) k is 1 for k = and it is less than 1 for k > 0, since n + a — 1 > 
n — 1 > n — k. For fixed k, our lower bound converges to 1 as n goes to infinity. Now, 
for all n > M we may write 

\n/2\ _ k M _ k L«/2J 



E 



P 



k\(a + k- 1 



■ o~n,k + 



k\{a + k- \)t 



k=0 y k=0 y k=M+l 

By our choice of M, and by \5 n ^\ < 1, the second sum on the right hand side has 
absolute value less than e/3, for all n > M. As n goes to infinity, each of 8 n ^, 
converges to 1, so there is an n such that for all n > n we have 



j On,M 



implying 



M 

E 

k=0 
\n/2\ 

E 

k=0 



k\(a + k-l) k 



k\{a + k-l) h 



M 



k=0 



k\(a + k - T) 



< 



>n,k 



M 

E 

k=0 



P 

kUa + k- 1) 



< 



2e 
3"' 



(2.2) 
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Combining (12. ip and (12.21) we obtain 



\n/2\ 

E 

k=0 



P 



k\{a + k-l) k 



P" 



k=0 



k\(a + k - l) k 



< e for all n > uq. 



□ 



Lemma 2.6. Let a be any positive real number and p be any real number. Then 



lim — 

n-voo (a + n — \) n 



L(n-1)/2J , i iM / , i ; \ 00 

(ra-fc-1)! /n + a- l-AA fc+1 = ^ 
Z^ ^1 i r Z^ 



m+l 



m!(cr + m) m +i 



k=0 v ' m=0 

Proof. The proof is very similar to the proof of Lemma 12. 5[ thus we omit the details. 
We only note that this time we have 



1 



(a + n- l) n 



^ (n-t-l)! /„ + „ - 1 - k\ t+1 _ l ^>< 

Z^ i.i I „ _ o h - 1 p z^ 



fc=0 



where 



n - 2k - 1 

(n-k- l) fc 



fc=0 



;i,A: 



(cr + n - l)fe 

Again \b~' nk \ < 1 for all n and and, for any fixed k, the limit of S' n k is 1 as n goes to 
infinity. □ 

Combining Lemmas 12.51 and 12.61 with Theorem 12.21 we obtain 



Pnd+r-l 



F r+1 (a)J2 



m=0 



ml(a + m — 1), 



+ (-l) d - r F d _ r . 1 (a)F d (a)f3 1 J2 



P 



m+l 



m=0 



m!(cr + m) m+1 



(2.3) 



and 



Qnd+r-l 



F d (a)"^(o- + n - l) n 



oo 

F r (a) £ 

m=0 



m!(cr + m — 1), 



-m+l 



m=0 



77l!((T + m) m+ i' 



(2.4) 

Taking the quotients of the right hand sides of (I2.3P and (12. 4p we obtain the following 
formula for £ = £(a, /3q, d, r): 



m=0 



m!(<T + m — l) r 



+ (-l) d -^F d _ r _ 1 (a)F d (a)/3i^ 



n m+l 



m=0 



m!(cr + m) m+1 



F r (a)J2 



P" 



+ {-l) d ~ r+l F d ^{a)F d {a)^ 



m=0 



ml (a + m — 1), 



OC 

E 

m=0 



/0 



m+l 



m!(cr + m) m+ i 



(2.5) 
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This is the formula that could have been discovered in a "parallel universe" where 
interest in Hurwitzian continued fractions were to arise a long time before defining 
Bessel functions. We may restate Eq. (12. 5p in terms of Bessel functions using the 
following two obvious statements. 

Lemma 2.7. For any positive p and a we have 

^m!(cr + m-l) m (M-i)/* aU ^ Q m\(a + m) m+1 ' p( CT ~ 2 )/2 

Here I v {z) is the modified Bessel function defined in U.6\) . 

Lemma 2.8. For any negative p and positive a we have 

P m J <T „ 1 (2 v ^)r(a) 
Z -m!(a + m-l) m (-p)C- 1 )/ 2 ana 

p m+1 J a {2^-p)T{a) 



J^ m \(a + m) m+1 (-p)(-2)/2 • 

Here J v {z) is the Bessel function given by (fi.Tp . 

Using Lemmas 12.81 and 12.71 above we may rephrase Eq. (12. 5ft as follows. 

Theorem 2.9. Let a, (3 , (3i andd be positive integers, and letr be nonnegative integers. 
Then the Hurwitzian continued fraction 



oo 
n=l 



£(ct:,/3o,/3i,cZ,r) = [ a, . . . , a , /3 , a, . . . , a, /3 + & ■ n 
is given by 

U ' J F r (a)/ CT _ 1 (2Vp) + (-irF d _ r («)/ CT (2 v ^) 

i/ d is odd, and it is given by 

F r+1 (a)J CT _ 1 (2 v ^) + (-l) r+1 F d _ r _ 1 («)J (T (2 v ^) 



F r (a) J (T _i(2V-p) + (-1) F d _ r \a)J a (2y/-p) 

if d is even. Here I v (z) and J v {z), respectively, denotes the is the modified respectively 
original Bessel function defined in hi. 6^1 and Hjl.7\) , respectively. 

Proof. We work out only the case of odd d in detail, the case of even d is completely 
analogous. Direct substitution of Lemma 12.71 into (12. 5p yields 

F r+M) (g _ 1)/2 + (-1) F d _ r _!(a)F d (a)^ )/2 

/3 ,/3i,d,r) =- 



pyu-L)/* p^ 

After multiplying the numerator and the denominator by p^ -1 )/ 2 we get 

F r+1 (a)/ CT „ 1 (2Vp) + (-l) d " r v ^F^ r _ 1 (a)F d (a)/3 1 / CT (2Vp) 



£(a,/?o,/3i,d,r) 



J P r (a)/ CT _ 1 (2Vp) + (-l)^' +1 VpF d _ r (a)F d (a)/3 1 J (T (2Vp) 
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The first equality in Theorem 12.91 follows from Jp~ = /3f i^(a) -1 and from the fact that 
(— l) d = (—1) in this case. A similar reasoning for even d yields 

, r R o , , _ F r+1 (a) J g _!(2y^) - (-l)^F rf _ r _ 1 (q)J CT (2 v ^) 
U a,p ,Pi,a,rj Fr{a) j a _ i{2V =- p) _ { -if-r^ Fd _ r{a) j ai 2^p) ' 

a final simplification may be made by observing that (— l) d = 1 in this case. □ 

Remark 2.10. Theorem 12.91 may also be derived from Lehmer's work [9] directly, using 
a few, easily verifiable facts about Fibonacci and Lucas polynomials. Again, we outline 
the proof for odd d only, the case of even d being completely analogous. Consider first 
the case when d — 1 and r = 0. In this case we get a = fio/fli, regardless of of a. 
Theorem 12.91 takes the form £(a,/3 ,/3i, 1,0) = J cr _i(2 v /p)// (T (2 v /p), which, by p = /3f 2 , 
is exactly Lehmer's formula (jl.5p . Consider next the case when d is an arbitrary odd 
integer and r = 0. In this case, we may use Lehmer's (9j Theorem 4] with all constants 
being equal to a. Using the fact that 

F d+ x(a) 

\a ,a\ - 



F d (a) ' 
d 

the fractions A/B, A'/B' and A" / B" appearing in |9] Theorem 4] are easily seen to 
correspond to ((/3 -a)F d (a) + F d+1 {a))/ F d {a), ((/3 -a)F d _i(a) + F d (a))/F d _ 1 (a;) and 
F d (a) / F d _i(a.) , respectively, in our notation. Lehmer's (bB+B'+B")/(aB) corresponds 
to our 

(3 F d (a) + 2F d _ 1 (a) = (ft - a)F d (a) + (aF d (a) + 2F d _ 1 (a)) 

PiF d (a) ^F d (a) ^ 

since aF d (a) + 2F d -\(a) = L d (a) holds for all d. Lehmer's [9J Theorem 4] gives 



F d (a) V v ' UmF d F d (a))^) 



which is the same as the formula implied by Theorem 12.91 after noting that we have 
p = (f3iF d (a))~ 2 . Finally, for arbitrary r we may substitute r] := £(a, /3q, fli, d, 0) into 



the formula 



and obtain 



F r+1 {a)r] + F r {a) 

\a,...,a,T]\- 



F r (a)7] + F r _i(a) 



W'h'h'W - Fia) ( i 7 Fa (a) , i^f^a - ia) 

P r\U) y Fd[a) y *d-l[.Ot) "I" / CT (2(/3 1 F d F d ( a ))-i) ) + *>-l W 

(Recall that a does not depend on r.) Equivalently, 

_ F r+1 (q)J CT _ 1 (2(/3 1 F d ( a ))- 1 ) + (F r (a)F d (a) - Fr+^F^ja)) J f7 (2(/3 1 F d F,(q))- 1 ) 
F r (a)/ CT _ 1 (2(/3 1 F (i (a))- 1 ) + (F r _i(a)F d (a) - F r (a)F d _!(a)) / (T (2(/3 1 F (i F d (a))- 1 ) ' 

Theorem 12.91 now follows after observing that F r (a)F d (a) — F r+ i(a)F d _i(a) and and 
F r _! (a)F d (a) — F r {a)F d _i{a), respectively, may be replaced by (—l) r+1 F d _ r _i(a) and 
(— l) r Fd-r(ot) respectively. 
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3. TWO USEFUL LEMMAS 

In this section we provide a combinatorial proof for two polynomial identities which 
seem to be interesting by their own right. They will play a crucial role in Section [4] where 
we prove Theorem 12.21 Our lemmas will be summation formulas for the polynomials 

1 ™ 

Rn(x,y) :=-V 

fc=0 

and 

*•<*• »> == JS S (*) (n -k - l) ( " " * )!( " " * " (3 ' 2) 
Lemma 3.1. The polynomials R n (x,y) satisfy 

t {=S*" (l " ») - s E (I) ("I' «*) (« " *)' V /- «« » * »■ 

m=0 v y fe=0 V 7 V 7 

Proof. We consider both sides of the equation as polynomials in the variable y with 
coefficients from the field Q(sc). Since a nonzero polynomial has only finitely many 
roots, it suffices to show that the two sides equal for any nonnegative integer value of 

V- 

For y G N, the left hand side may then be rewritten as 

^t( B :J(-)-i:(:)(;;fl(".-*)>v. 

m=0 v 7 fc=0 V 7 Vy 7 

This is 1/n! times the total weight of all quadruplets (^1,^2,71,72) subject to the 
following conditions: 

(1) 7Ti is a permutation of {1, . . . , n}, 712 is a permutation of {1, . . . ,n + y}, the 
functions 71 : {1, . . . , n) — > {0, 1, 2} and 72 : {1, . . . , n + y} — > {0, 1, 2} are 
colorings; 

(2) for % G {1, 2} an element j satisfying = 1 must be a fixed point of 7Tj; 

(3) for any j, 71 (j) = 2 is equivalent to 72(7) = 2 and 71 (j) = 72O) = 2 implies 
that j is a common fixed point of 7Ti and 7r 2 ; 

(4) the colorings 71 and 72 satisfy \{j e {1, . . . , n + y} : 7 2 (j) > 0}| = \{j G 
{l,...,n} : 71(7) > 0)1+2/; ' 

We define the weight of (771,772,71,72) as 

a .|{je{i J ...,n}:Tn(j)=i}| . r_ x \\{je{i,~,n}--7iti)=2}\_ 

Indeed, after setting n-mas the number of elements j satisfying 71 (j) = 72 (j) = 2, 
there are ( " ) ways to select them. By rule ([3]) these are common fixed points of 7Ti 
and 7r 2 (and thus elements of {1, ... , n}). Next we set k as the number of elements j G 
{1, . . . ,n} satisfying 71 (j) = 1, and select these elements, in (™) ways. The remaining 
elements of {l,...,n} satisfy 71 (j) = 0. At this point we have \{j G {1, . . . ,n} : 
7i(j) > 0}| = n — m + k. Thus, by rule (J1J), the set {j G {1, . . . , n + y} : 7 2 (j) > 0} 
must have n — m + k + y elements. Exactly n-mof these elements satisfy 72 (j) = 2, 
therefore the number of elements j G {1, . . . , n + y} satisfying 7 2 (j) = 1 must be y + k. 



(:)(:-:>-*>** 
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There are (™+^ ) ways to select them. So far we have selected n — m + k fixed points of 
7Ti and n — m + y + k fixed points of 7T2. We can complete prescribing the action of 7r± 
and 7r 2 in (m — k)\ 2 ways. 

The same total weight may also be found by fixing the pair (jri, 7r 2 ) first, and summing 
over all allowable pairs of colorings (71,72)- We call the pair (71,72) allowable, if the 
quadruplet (7Ti, 7r 2 , 71, 72) satisfies the conditions listed above. Suppose j is a common 
fixed point of 7Ti and i\2- Observe that the contribution of all allowable pairs (71,72) 
satisfying 71 (jo) = 72 (jo) = 2 cancels the contribution of all allowable pairs satisfying 
7i(jo) = 72(jo) = 1- Indeed, let (71,72) an allowable pair satisfying 71 (j ) = 7 2 (j ) > 0. 
Then the pair (7^, 7 2 ) defined by 

iti) - 

is also allowable and also satisfies 7i(jo) = I2U0) > 0. (Here j G {1, . . . ,n} for 71 and 
7i and j G {1, . . . , n + y} for 72 and j' 2 .) The map (71, 72) !->■ (7^ 7 2 ) is an involution 
that matches canceling terms: the only difference between the respective contribution 
is a factor of x or — x associated to jo- Therefore we may restrict our attention to the 
total weight of quadruplets (^1,^2,71,72) satisfying the following additional criteria: 

(5) no j satisfies 71 (j) = 72 (j) = 2, in particular, the colorings 71 and 72 map into 
the set {0, 1}; 

(6) no j satisfies 71 (j) = 72 (j) = 1, in other words, the sets {j G {1, . . . , n} : 
7i(j) = 1} and {j G {1, . . . , n + y} : j 2 (j) = 1} are disjoint. 

When computing the total weight of such quadruplets, we may first select k as the 
number of elements j G {1, . . . , n} satisfying 71 (j) = 1 and select them in (™) ways. 
By rule (jl]) there must be y + k elements j G {1, . . . , n + y} satisfying J2U) — 1 an d, 
by rule (6), this set is disjoint of the previously selected fc-element set. Thus there are 
("y+fc ) ways to select them. So far we have selected k fixed points of 7Ti and y + k fixed 
points of 7r 2 . There are (n — k)\ 2 ways to complete prescribing tti and 7r 2 , the weight 
of the quadruplet is x k . We obtain a total contribution of Y^l=o Cy+k K n ~ k)\ 2 x k 
which is exactly n\ times the right and side. □ 

Lemma 3.2. The polynomials S n (x,y) satisfy 

i ■ h l xT (;) (;_T-x) <» - *>■<- - * - ^ 

for all n > 0. 

Proof. A proof may be obtained by performing slight modifications to the proof of 
Lemma [3. 1\ which we outline below. The left hand side equals 



s £ L " J ( -^™^ (T) LTkl 1 ) (m - mm - k - i)]x 



m=0 x / fc=0 



k 



which may be considered as x/n\ times the total weight of quadruplets (^1,^2,71,72) 
where the only change to the definition is that, instead of (4), now we require 
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(4') the colorings 71 and 72 satisfy \{j G {1, . . . , n + y} : 72(7) > 0}| = \{j G 
{l,...,n} : 7 i(j) >0}|+y + l. 
Thus we will have to select y + k + 1 elements j (instead of y + k) satisfying vr 2 (j) = 1, 
in ( ^t^J • In the last stage, we will have selected n — m + y + k + 1 fixed points of 7T2, 
thus we will only have (m — k — 1)! ways to complete the selection of 7T2. 

The same involution as before shows that we may again restrict our attention to those 
quadruplets which satisfy the additional conditions (5) and (6). Let us set k again as 
the number of elements j G {1, . . . ,n} satisfying 71 (j) = 1. The only adjustment we 
need to make to the reasoning is to observe that now we need to have y + k + 1 elements 
satisfying |7 2 (j) = 1| which may be selected in ways, instead of ways. 

Finally, we may complete the selection of 7T2 in (n — k — l)\ ways, instead of {n — k)\. The 
resulting total weight is exactly x/n\ times the right hand side of our stated equality. □ 

For the sake of use in Section 0] let us note that R n (x,y) and S n (x,y) may also be 
written in the following shorter form, using falling factorials: 

fc=0 ' k=0 

A similar simplification yields that Lemma 13.11 is equivalent to 

m=0 v y fe=0 v 7 

and that Lemma [3.21 has the compact form 

± V ) = 1< "E " ( " + i * ) for all n > 0. (3.5) 

^^m — ml! k\ \n — 2k — 1 J 

m=0 v y fc=0 v 7 

4. Calculating the convergents 

In this section we calculate (nd + r — l)st convergent of £(o:, flo, fli, d, r) directly, 
from the Euler-Mindig formulas (11. 3p . We begin by observing that, for an arbitrary 
continued fraction [do, ai, . . .] and an arbitrary positive integer a, the numerator p n in 
(11. 3p may be rewritten as 

Pn = ^(fli -« + «)= X] JJ(^i - «)a'^ T| - 

5C e {0,...,n} ieS 5C e {0,...,n} TCS ieT 

Changing the order of summation gives 

TC{0,...,n} iST TCSC e {Q,...,n] 

For £(a, (3q, 0i, d, r) we have cii = a unless i is congruent to r modulo d. Thus, to 
compute p n d+r-i using (14.11) . we only need to sum over subsets T whose elements are 
all congruent to r modulo d. This observation yields the following recurrence: 

n-l 

Pnd+r-i = F nd+r+1 (a) + J^pw+r-i • (A> + Pi ■ k - a) ■ F {n _ k)d {a). (4.2) 

fc=0 
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By Eq. (11. lip , the term F nd+r+ i(a) above is the contribution of T = 0, whereas the 
term pkd+r-i • (A) + fti-k — a) • F( n _ k ) d+1 (a) is the total contribution of all sets T whose 
largest element is kd + r. Substituting n = into (14. 2 p yields the initial condition 

Pr-i = F r+1 (a) 

which is obviously true for r > 0, and it is also valid when r = after setting p_i = 1, 
as usual. Consider the formal Laurent series 

oo 

y(t) := J^Pnd+r-i ■ t^ n ~ a e Q((t)). 

n=0 

For y(t), Eq. §4~2§ yields 

oo oo 

y(t) := F nd+r+1 (a)t^ n ■ t^~ a + 1 ■ ^ F nd {a)t^ n ■ y'(t), (4.3) 

n=0 n=l 

where is the formal derivative of y(t) with respect to t. To write (14. 3 p in a more 
explicit form, observe that, by Eq. ( I1.12p . we have 

00 00 n nd+r+l „nd+r+l i / _r+l _r+l \ 

Using the fact that p\p2 = —1, the above equation may be rewritten as 

^ WJ v^n i-( P f + ^)t + (-i)^ 

1 ^+1 _ p r+l + ( _l)r+l ( ^-r-l _ p ^-l )t 



v^+4 l-(^ + p|)t+(-l)^ 2 

By Eqs. (011 and 1Q3J) we obtain 



= i_ Ld(Q)t+( _i)«i t 2 • ( 4 - 4 ) 
Note that substituting r = d — 1 in Eq. ( 14.41) yields 

n=l n=0 

since .Fo(cc) = 0. Using these last two equations, we may rewrite (14.31) as 

F T+1 {a) + (-^) r+l F d -r-i{a)^ 0n _ a F d (a) ■ 

g[ ' 1-L d (a)tfr + (-1)^ + l-Li(a)tft + (-l)^ ^ U ' 

Rearranging to express the derivative of y(t) yields 

y(t) ■ (r*- 1 - L d (a)t-i + (-l)^- 1 ) 



(F r+1 (a) + (-l)^.^^)^) ■ t*°- 



(4.5) 
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Inspired by the way we solve ordinary differential equations in analysis, we will guess 
the solution of ( 14. 5 p by "solving" first the corresponding homogeneous equation and 
then replace the arbitrary constant by a formal Laurent series. The next few lines will 
not make sense, they just indicate how one may come up with a good guess for y(t). 
A reader who does not like "obscure reasoning," should skip ahead to (14.61) and accept 
that there is no "rational explanation" as to why introducing the formal Laurent series 
z(t) is a good idea. 

The homogeneous equation 

, m _ ^(t)-(r^- 1 -L d ( Q )t- 1 + (-i)^- 1 ) 
VH[t) ~ TM) 

"may be rewritten as" 

d ... r*- 1 - u^t- 1 + (-l)^- 1 

\n(y H (t)) 



z(t):= y(t) ■ 



exp 



dt yany,J F d (a) 

"yielding" 

V U\ _ n . f -L d (a)/F d (a) ( + (-l)^ 1 

where C is an arbitrary constant. This "solution" to the "homogeneous equation" 
suggests looking for a solution to Eq. (14. 5 j) of the form 

y (, )=;(l ),^ Wtip ( -';;^j» 

Equivalently, we would want to set 

-Pi + (_i)^i^i 
Ji ■ F d {a) 

Alas, the resulting formal expression would contain arbitrary large positive and as well 
as arbitrary small negative powers of t, it does not resemble a Laurent series at all. 
Hoping that a slight change would not upset our calculations irreparably, we define z(t) 
by setting 

z(t):= y{t) . t , m («) 

Note that z(t) is a formal Laurent series in the variable tV-^w ; an infinite formal sum 
of the form 

oo 

z(t) = J2s n t^ n+f}o ~ a+Ld[a)/Fdia) e Q{{t 1/Fd(a) )). 

n=0 

Taking the derivative on both sides of ( 14. 6 j) yields 

lm) f{-l) d -H Pl 
z '(t) = y'(t) ■ t*SW exp 



Pi ■ F d (a) 

+y(t) ■ exp ' 



f3 1 -F d (a)J V F d {a 
exp \ . r • y\t) + y(t) 
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(Note that L d (a) / F d (a) is always positive.) After substituting the value of y'(t) from 
( 14. 5 p and simplifying we obtain 

z\t) = exp V(t) ~ (*h-i(«) + (-lr^^a)**)**- 



0i ■ F d (a) J F d (a)t^ 
By (I4.6p . the last equation is equivalent to 

z(t) 



— t F d(«) exp 1 



Comparing the coefficients of t^ i n+/3 ° a + L d.w/ F dW i on D0 th sides yields 
/?i ■ n + A) - a + 



(-1) 

L d (a)\ s n+l F r+ i(a) \h-F d (a) 



F d (a)J F d (a) F d (a) (n + 1)! 

+ (-l) r F d _ r _i(a) ^/?i-F d (a) 



F d (a) n! 

Note that the factor ((3i-n + (3 — a + L d (a)/ F d (a)) in the last equation equals (3i-(n+a), 
where a is the magic sum, and that the magic quotient p appears twice on the right 
hand side. Since, by Corollary 12.41 (c + n )n+i is not zero, we may divide both sides by 
F d (a) n (3[ l+1 (a + n) n+ i, and obtain the following recurrence for s n := s n / (F d (a) n [3i(a + 
n - l) n ): 

s n+1 = s n + F r+ i(a) • 



(n + l)!(a + n) n+ i ^(a)^ n!(a + n) n+ i' 

Considering the fact that So = s = F r+ i(a) and that (— l) r F d _ r _i(a) / (F d (a) f3i) 
(— l) d ~ l ~ r F d _ r _i(a)F d (a)f3i ■ p, the last recurrence implies 

n k n— 1 

*• = ■ E M(g A-i )t + (-i)~*~.<°w«)A • E 

/c — — 

Multiplying both sides by (er + n) n+ i yields 

Sn „ , N p fc (^ + n) n _ k 

wff - " +l(a) '^ « 

^d-r-zP /^n/ Afl ^p^ +i (a + ™)n-fc-i 



+ (-l) d - r F d _ r _ 1 (a)F d (a)/3 1 -^ 



fc=0 



A'! 



By the formulas given in Eq. 03. 3p . we may rewrite the previous equation in terms of 
the polynomials R n (x,y) and S n (x,y) as follows: 

^— = F r+1 (a)R n (p, a - 1) + (-^-^^(a)^)/?! • S n (p, a - 1) . (4.7) 



F d (a) n P[ 
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By (14.61) we have 



y {t) = z(t) ■ r^W/^W exp 
Comparing the coefficients of tpo+pi-n-a yi e ici s 



Pnd+r-1 — ^ ] 
m=0 

hence we have 



(-i) a 



(n — my. 



Pnd +r-l _ V- S m ^F d (q)2/3f ^ _ ^ S m p 



jn—m 



F d (a) n P? ^ F d (a) m ^ (n-m)! ^ F d (a)^/3f (n-m)!' 

Substituting Eq. (14. 7p into this last equation and using Eqs. (13 .4p and ( 13. 5 p yields the 
formula for p nd+r _i in Theorem 12.21 

For positive r, the formula for q nd+r _i stated in Theorem 12 .21 is an easy consequence of 
the formula for p n d+ r -i- Indeed, by Lemma \L2\ the denominator q n d+ r -i is the same as 
the numerator p nd+r _ 2 associated to £(a, ft, ft, d,r — i), thus we only need to replace r 
by r — 1 in the formula stated for p nd+r _ 2 . It only remains to show the following lemma. 

Lemma 4.1. The equation stated for q nd+r _i m Theorem \2.S\ remains valid when we 
substitute r = 0. 

Proof. Since Fo(a) = 0, substituting r = in Theorem 12.21 gives 



L(n-1)/2J , 7 i \ l / i ;\ 

1} Fd{a) 2. jfcj { n -2k-l ) 



p h+1 . 



in — k — 1)! fn + a — 1 — k\ 

fc=0 



fc=0 

Here we may replace ( — l) d+1 F d (a) 2 ■ ft by p^ft" . After rearranging we obtain 

L(n-1)/2J 

^n — 

~ fc! V " - 2fc - 1 
We need to show the validity of this equation. Observe that 

£(a,ft,ft,d,0) = [ft, a, . . ■ , a, ft + ft, . . .] = ftp - a + [ a, . . . , a , ft + ft, ■ ■ .] 

= ft - a + £(a:,ft + ft,ft,d,d) 

Thus the denominator g n( f_i associated to ft, ft, d, 0) is the same as the same as 
the denominator q(n — l)d — 1 associated to £(a, ft + ft, ft, cZ, d). We may apply the 
already shown part of Theorem 12.21 Since the magic quotient p depends only on a, ft 
and d, it is the same for ft, ft, d, 0) and for £(a, ft + ft, ft, d, <f). For the magic 
sums we get 

( a , a o j j\ ft + ft ~ a . L d (a) 

a{a, ft + ft, ft, a, a) = h - , . = <r(a, ft, ft, d, 0) + 1. 

Pi Pli'dK.Oi) 
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Therefore we may obtain an equation for the q n d-i associated to £(a, /?o, d, 0) by 
replacing n with n — l,r with d and cr with a + 1 in the formula for q n d-\ in Theorem 1 2. 2 1 
Since Fo(a) = 0, the second sum vanishes and we get 



qw+r-i(a, A),/#i,d,0) p / L(r L^ /2j (n - 1 - k)\ (n + a - 1 - k 

~ Fd{a) 2. jfei ^ n _i 



which is obviously equivalent to (14. 8p . □ 

5. Special cases leading to elementary expressions 

In this section we describe all instances of Theorem 12.91 for which the magic number 
a is the half of an odd integer, forcing all (modified) Bessel functions in the statement 
to be known elementary functions. There is not much to say about the case d — 1: as 
seen in Remark 12. 10} we have a — /?o/A thus a depends only on (3 and fli in a very 
simple fashion. In the case when d > 2 we give a similarly simple description. 

Theorem 5.1. If d > 2 then o is the half of an odd integer if and only if one of the 
following conditions holds: 

(1) d = 3, a = 1 and (/3q + is the half of an odd integer; 

(2) d = 2, a = 1 and (/?o + 2)//3i is the half of an odd integer; 

(3) d = 2, a = 2 and ((3q + t/ie /ia// o/ an odd integer; 

(4) d = 2, a = 4 and (2/3q + l)//3\ is an integer. 

Proof. First we show that even assuming that a is the half of an (even or odd) integer 
implies that d is at most 3. The fact that a belongs to (1/2) • Z implies the same for 
Pier + a — /3q = Ld(ot)/ Fdipt). We may subtract any integer multiple of Fd(a) from 
Ldipi) and still have element of (1/2) • Z. Let us select r^(a) := Ld(a) — aFd(a) and 
consider the fraction ra(a) / Fd(a) G (1/2) ■ Z. We claim that 

< ^4^1 < 1 holds for a > 2 and d > 3. (5.1) 

Indeed, for d = 3, we have r 3 (a) = L 3 (a) — aF 3 (a) = (a 3 + 3a) — a(a 2 + 1) = 2a and 
F 3 (a) = a 2 + 1. Both r 3 (a) and -F 3 (a) are positive and r 3 (a) < -F 3 (ct) follows from 

-FMaO — r 3 (a) = (a — l) 2 > for a > 2. 

For d = 4 we have r 4 (a) = L 4 (a) — «F 4 (a) = a 4 + 4a 2 + 2 — a(a 3 + 2a) = 2a 2 + 2 and 
-F4(«) = a 3 + 2a. Both r±(a) and ^4(0:) are positive and r^{a) < F^{a) follows from 

F 4 (a) - r 4 (a) = a 3 - 2a 2 + 2a - 2 = a 2 {a 1 - 1) + 2(a - 1) > for a > 2. 

For larger values of d, we may show that 

< r<i(a) < Fd(a) holds when d > 3 and a > 2 

by induction on d, using the fact that the statement is valid for d G {3, 4} and that 
rd(a) satisfies the same recurrence as Fd(a), allowing to express the inequality for the 
next value of d as a positive combination of the inequalities for the current and the 
previous values of d. This concludes the proof of (15. ip . 

As a consequence of the inequality (15. ip . whenever d > 3 and a > 2, the only way 
for rd(a) / F d (a) to be an integer is to have r d (a)/ F d (a) = 1/2. In other words, in this 
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case, we must have 2r d (a) = F d (a) which is equivalent to 2L d (a) = (2a + l)F d (a). To 
show that a > 2 and d > 3 can not hold simultaneously, it suffices to show that 2L d (a) 
can never be equal to (2a + l)F d (a). To do so, first we observe that 

2L d (a) > (2a + l)F d (a) holds for a G {2, 3} and d > 3. 

Indeed, for a = 2 we have 28 = 2L 3 (2) > (2 ■ 2 + 1)F 3 (2) = 25 and 68 = 2L 4 (2) > 
(2 • 2 + 1)^4(2) = 60 and we may prove the same inequality for higher values of d by 
induction. Similarly, for a = 3 we have 72 = 2L^(3) > (2 ■ 3 + 1)^3(3) = 70 and 
238 = 2/74(3) > (2 ■ 3 + l)i<4(3) = 231 and we may proceed again by induction on d. 
Finally, to exclude a > 4, we will show 

2L d (a) < (2a + l)F d (a) holds for a > 4 and d > 3. 

For d = 3 we have 

(2a + l)F 3 (a) - 2L 3 (a) = a 2 - 4a + 1 = a(a - 4) + 1 > 0, 

and for d = 4 we have 

(2a + l)F 4 (a) - 2L 4 (a) = a 3 - Aa 2 + 2a - 4 = a 2 (a - 4) + 2a(a - 2) > 0, 

and again we may proceed by induction on d. 

We obtained that, for d > 3, a G (1/2) ■ Z is only possible if a = 1. In that case, 
for d > 4 we have L d (l) — 2F d (l) = F d _ 3 (l) (this may be shown by induction). Clearly 
L d (l) / F d (l) is the half of an integer, if and only if the same holds for F d _s(l) / F d (l) . 
Now we may exclude the case d > 4 completely, after observing that 

Fj qH") 1 

< d ~ 3 \ J < - holds for d > 4. 
^(1) 2 

Indeed, the above inequality is equivalent to < 2F d ^ 3 (l) < F d (l) which may be easily 
shown by induction. 

We have shown that a can only be the half of an integer if d = 2 or d = 3. In the 
case, when d = 3, we have also shown that only a = 1 is possible, and we get 

_ A)-i £ 3 (i) A)-i _±_ Po + l 

° ft ^F 3 (l) ft 2ft ft ' 

Consider finally the case when d = 2. As before, a G (1/2)Z implies that L d (2)/F d (2) = 
(a 2 + 2)/a G (1/2)Z. This implies that a must be a divisor of 4, that is, an element of 
{1,2,4}. We have 

ft — ol L 2 (a) aft — a 2 a 2 + 2 aft + 2 

ft ft-^McO ft« fta; fta 

Therefore, for a = 1 we get cr = (ft + 2)/ft, for a = 2 we get a = (ft + l)/ft and for 
a = 4 we get a = (2ft + l)/(2ft). □ 

A nice example of the case when d = 3 and a = 1 in Theorem 15.11 above is the case 
when a = 1, ft = 3m — 1, ft = 2m, = 3, r = 2 for some m > 0. In this case we get 
a = 3/2 and p = l/(16m 2 ). Theorem 12.91 gives 

2/ 1/2 (l/(2m)) 



f(l,3m- 1,2m, 3, 2) 



J 1/2 (l/(2m))+J 3/2 (l/(2m)) 
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Using (11. 8ft and (ll.lOjl we may rewrite the preceding equation as 

2sinh(l/(2m)) 



[1,1,3m- 1 + 2mn\™ 



cosh(l/(2m)) - (2m - 1) sinh(l /(2m)) 
For m = 1 we obtain 



[1,1,2,1,1,4, 



e l/2 _ e -l/2 



e l/2 

A similarly nice example for the case when d = 2 and a = 1 in Theorem 15.11 above is 



f (1, 3m - 2, 2m, 2, 1) = [1, 3m - 2 + 2mn]~ 
for some m > 0. In this example a = 3/2 and p = (— l)/(4m 2 ) hold. Theorem 12.91 gives 

Ji /2 (2v^p) J 1/2 (l/m) 



J 1/2 (2 A / = p) - J 3 / 2 (2v /3 p) Ji/ 2 (l/m) - J 3 /2(l/m)' 

Using ( II. 9p and (jl.lOp we may rewrite £ above as 

sin(l/m) 
cos(l/m) — (m — 1) sin(l/m) 

Substituting m = 1 yields tan(l) = [1, 1, 1, 3, 1, 5, 1,7,...]. 

We conclude this section with an example which is not likely to be found in the 
literature, due to its "sheer ugliness". Let us set a = 4, (3q = 7m + 3, f3\ = 2m + 1, 
d = 2 and r = 1, where m is any nonnegative integer. For this example we have a = 7/2 
and p = (— l)/16(2m + l) 2 , yielding l^^p = l/(4m + 2). Using the fact that 

J 5/2(z) = ((J^-l) sin(z) - ^cos(z) 

and 

■m*) = V^((S"i) Mz) - (5 " cos(z) 

(see PHI List of formula?: 46, 47]), Theorem 12.91 gives 

4 ((12(2m + l) 2 - 1) sin (^) - 6(2m + 1) cos (^)) 



(240m 2 + 228m + 53) cos (j^) - (960m 3 + 1392m 2 + 648m + 97) sin (4^) ' 



The above formula was calculated from Theorem 12 . 9 1 with the help of Maple. The same 
program was used to double-check its correctness for selected values of m. For example, 
for m = we obtain 

4 (11 sin (i) -6cos(i)) . 

-i TTT^ ^ = 4, 3, 4, 4, 4, 5, 4, 6, 4, 7, 4, . . . . 

53 cos (|) -97sin(±) L '''''''''' ' 
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